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Abstract 



Transfection of an anionic polynucleotide through a negatively charged 
membrane is an important problem in genetic engineering. The direct associ- 
ation of cationic surfactant to DNA decreases the effective negative charge of 
the nucleic acid, allowing the DNA-surfactant complex to approach a nega- 
tively charged membrane. The paper develops a theory for solutions composed 
of polyelectrolyte, salt, and ionic surfactant. The theoretical predictions are 
compared with the experimental measurements. 

PACS.05.70.Ce - Thermodynamic functions and equations of state 
PACS.61.20.Qg - Structure of associated liquids: electrolytes, molten salts, etc. 
PACS.61.25.Hq - Macromolecular and polymer solutions; polymer melts; swelling 
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Gene therapy has increasingly captured pubhc attention after the first gene transfer study 
in humans was completed in 1995 The procedure delivers a functional polynucleotide 

sequence into the cells of an organism affected by a genetic disorder. The gene delivery 
system that has been adopted in over 90% of the clinical trials to date is in the form of ge- 
netically engineered non- replicating retroviral or adenoviral vectors 0. Unfortunately, the 
adverse response of the immune system has hindered application of virus based gene therapy. 
New strategies are now being explored [^|T9(]. One of the approaches, pioneered by Feigner 
and Ringold , relies on association between the anionic nucleic acid and cationic lipid lipo- 
somes. The process of association neutralizes the excess negative charge of a polynucleotide, 
allowing the DNA-lipid complex to approach a negatively charged phospholipid membrane. 
Unfortunately, the cationic lipids and surfactants are toxic to an organism. A question that 
we will try to answer in this letter is: What is the minimum amount of cationic surfactant 
or lipid that is necessary to form a complex and how does this amount depends on various 
properties of a system? 

We study a solution consisting of DNA segments of density Pdna, surfactants of density 
Psurf, and salt molecules of density psait The solvent is idealized as a uniform medium 



of dielectric constant D. Since the DNA molecule has a large intrinsic rigidity, we model it 
as a cylinder of fixed length and diameter. When in solution, the Z phosphate groups of 
the DNA strand become ionized, resulting in a net molecular charge, —Zq. An equivalent 
number of counterions of density ZpoNA are released into solution preserving the overall 
charge neutrality. Similarly, the cationic surfactant molecule in aqueous solution becomes 
ionized, producing a free negative ion and a flexible chain consisting of one positively charged 
hydrophilic head group and a neutral hydrophobic tail. The ions of salt, the counterions, 
and the negative ions dissociated from the surfactant are modeled as hard spheres with point 
charge located at the center. For simplicity, we shall call the negative ions, "colons", and 
the positive ions, "counterions" — independent of the species from which they are derived 
(see Figure 1). 

The strong electrostatic attraction between the counterions, cationic surfactant, and the 



DNA favors formation of clusters consisting of one DNA molecule and ricount associated 
counterions, and rigurf associated surfactants. The process of association neutralizes rigurf + 
'T'count phosphate groups of a DNA molecule, decreasing the net charge of a complex to, 

^complex 

— —{Z — risurf — JT-count)? (Figurc 1). Our task is to determine the values of ricount 
and rigurf which are thermodynamically favored, i.e. which minimize the overall Helmholtz 
free energy of solution. 

For a dilute suspension, the main contributions to the free energy can be subdivided into 
three parts: the energy that it takes to construct an isolated complex, ^association; an energy 
that it takes to solvate this complex in the ionic sea, ^solvation; and the entropic energy of 
mixing, F^ixing. 

To calculate the free energy of an isolated cluster, we use the following simplified model 
of a complex. Each monomer of a polyion is treated as free or occupied by a counterion or 
a surfactant (Figure 2). We associate with each monomer i two occupation variables a{i) 
and T(i), such that a{t) = 1 if the site is occupied by a condensed counterion, and a{i) — 
otherwise. The occupation number for surfactants, T(i), behaves in a similar way. The 
free energy can now be calculated as a logarithm of the Boltzmann sum over all possible 
configurations of condensed counterions and surfactants along the polyion, 

/^-^association 1^ ^ ] € ^ . ( 1 ) 

The energy of a given configuration v can be subdivided into two parts, E^, = Ei + E2, 
where the electrostatic contribution is, 

'~^h D\r{i)-r{j)\ ■ 

The energy E2 arises from hydrophobicity of surfactant molecules. Clearly, when two ad- 
jacent sites are occupied by surfactants, the net exposure of hydrocarbon tails to water is 
reduced. We capture this effect by introducing an additional contribution to the overall 
energy of interaction, E2, given by 

^2 = -f E r{i)T{j) , (3) 

<i¥^j> 
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where the sum runs over the nearest neighbors. The parameter x, related to the decrease 
in overall energy due to the agglomeration of surfactants, is obtained from an independent 
experimental measurement of the energy that is required to move one surfactant molecule 



from a monolayer to bulk [BO] 



The exact solution of even this one-dimensional "sub-problem" is rather difficult due 
to the long-ranged character of the Coulomb force. To proceed we could use a mean-field 
approximation, but while the mean-field theory works very well for long-range potentials, for 
one-dimensional systems with short-range forces, it can lead to unphysical instabilities ^1 



In order to avoid this difficulty, we treat the long-range electrostatic part of the association 
free energy using a mean-field approximation, while performing an exact calculation for the 
short-range hydrophobic interaction. 

Once a cluster, constructed in isolation, is introduced into solution, it gains an additional 
energy due to electrostatic interactions with the other entities. The free energy gained in 
the process of solvation can be obtained using the Debye-Hiickel [ pT| , |22| theory. Let us fix 
the position of one cluster and ask what is the electrostatic potential $ that this cluster feels 
as a result of the presence of all the other clusters, surfactants, counterions, and colons. To 
answer this question, consider the Poisson equation, 

= -^P, . (4) 

To make the problem well posed, this equation must be provided with a closure which would 
relate the electrostatic potential, $, to the net charge density, pg. A simple closure motivated 
by ideas derived from the Debye-Hiickel theory is to suppose that the free (unassociated) 
surfactants, counterions, and colons are distributed around the complex in accordance with 
the Boltzmann distribution, with other clusters providing a neutralizing background, 

Pq = gcomploxPDNA + gPcountC"^'^* - gPcoione^^'^* + gPsurfC"^^*, (5) 

where /? = l/^ksT). 

Inserting (5) into (4) we obtain the Poisson-Boltzmann equation, which after linerization 



reduces to the familiar Helmholtz form. The hnearization is justified since all the non- 
linearities are effectively included in the renormalization of DNA charge by the formation of 
clusters p3|-p5| . The Helmholtz equation can be solved analytically, yielding the electrostatic 
potential of a complex. The electrostatic free energy of solvation is obtained through the 
usual Debye charging process |^,^ . 



The free energy due to mixing of various species is a sum of individual entropic contri- 
butions, 

Fmixing Fcounterion ~l~ Fcoion ~l~ Fsurfactant ~l~ Fcomplex- (6) 

The structure of each one of these terms is similar to that of an ideal gas and can be 



calculated using the Flory's theory of polymer melts p6|. For example, in the case of 
colons, the reduced free energy density is PFcoion/V = Pcoionln(0coion/C) - Pcoion, where 
0coion is the volume fraction occupied by the negative microions, while C is a factor that 
takes into account the internal structure of each specie. For structureless particles such 
as colons and counterions, C = 1- For ^ flexible linear surfactant chain, ( is the number of 
monomers comprising a molecule . For a complex made of a DNA segment and condensed 
surfactants and counterions, ( is related to the number of different configurations which can 
arise when ricount counterions and Ugurf surfactant molecules associate to a DNA molecule. 
Minimizing the total free energy, 

F -^mixing ~l~ -^association ~l~ -^solvation; ("i^) 

with respect to the number of associated counterions and surfactants, 

OF OF 

0, (8) 



we find the thermodynamically preferred values for the number of condensed particles. We 
define a "surfoplex" to be a complex in which almost all of the DNA's phosphate groups 
are neutralized by the associated surfactant molecules. As mentioned in the introduction, 
we are interested in the minimum amount of cationic surfactant needed to transform naked 



DNA into surfoplexes. To this effect, we study the dependence of the number of condensed 
surfactant molecules, ngurf, on the bulk concentration of surfactant. Figure 3 demonstrates 
the location of the cooperative binding transition associated with the formation of surfo- 
plexes. The transition is a result of competition between the entropic, the electrostatic, and 
the hydrophobic interactions. Clearly, for small concentrations of surfactant, binding to the 
polyions is not thermodynamically favored, since the system can lower its free energy by 
keeping the surfactant in the bulk, thus gaining entropy. As the density of surfactant in- 
creases, the gain in electrostatic energy, due to binding, outweighs the loss of entropy due to 
confinement which, in any case, is largely compensated by the release of bound counterions 



IS 



27| , |28[| . The cooperativety predicted by our theory and observed in experiments p9| , |30 
due to the fact that once the first surfactant is bound, the binding of additional surfactants 
is strongly favored by the decrease in hydrophobic energy of the exposed hydrocarbon tails. 
The high degree of cooperativety allows us to clearly define how much surfactant is necessary 
to form a surfoplex. To check the predictions of our theory, we compare it to the recent 



experimental measurements of binding isotherm in a solution of DNA, dodecyltrimethy- 
lammonium bromide, and salt, (Figure 3). The agreement is encouraging specially since the 
theory does not have any fitting parameters. 
ACKNOWLEDGMENTS 

We would like to acknowledge helpful conversations with Profs. K.A. Dawson, Michael 
E. Fisher, and Prof. H. E. Stanley. This work was supported in part by CNPq - Conselho 
Nacional de Desenvolvimento Cientffico e Tecnologico and FINEP - Financiadora de Estudos 
e Projetos, Brazil. 



6 



REFERENCES 
[1] R. M. Blaese et al, Science 270, 475-480 (1995). 

[2] M. J. Hope, B. Mui, S. Ansell, and Q. F. Ahkong, Mol. Membr. Biol. 15, 1-14 (1998). 
[3] W. F. Anderson, Nature 392, 25-30 (1998). 

[4] P. L. Feigner, T. R. Gadek, M. Holm, M., Roman, M. Wentz, J. R Northrop, M. 
Ringold, and H. Danielsen, Proc. Natl. Acad. Sci. U.S.A. 84, 7413-7417 (1987). 

[5] P. L. Feigner, and G. M. Ringold, Nature 337, 387-388 (1989). 

[6] T. Friedmann, Sci. Am. 276, 80-85 (1997). 

[7] P. L. Feigner, Sci. Am. 276, 86-90 (1997). 

[8] R. Langer, Science 249, 1527-1533 (1990). 

[9] I. M. Verma, and N. Somia, Nature 389, 239-242 (1997). 

[10] S. T. Crooke, and C. F. Bennett, Annual Review of Pharmacology and Toxicology 36, 
107-129 (1996). 

[11] C. P. Hodgson, Biotechnology 13, 222-225 (1995). 

[12] S. Akhtar, and S. Agrawal, Trends in Pharmacological Sciences 18, 12-18 (1997). 

[13] C. F. Bennett, D. Kornbrust, S. Henry, K. Stecker, R. Howard, S. Cooper, S. Dutson, 
W. Hall, W. and H. I. Jacoby, Journal of Pharmacology and Experimental Therapeutics 
280, 988-1000 (1997). 

[14] S. M. Sullivan, Journal of Investigate Dermatology 103, 85S- 89S (1994). 

[15] R. L. Juliano, S. Akhtar, Antisense Research and Development 2, 165-176 (1992). 

[16] X. Gao, and L. Huang, Gene Therapy 2, 710-722 (1995). 

[17] J. G. Lewis, K. Y. Lin, A. Kothavale, W./M. Flanagan, M. D. Matteucci, R. B. De- 

7 



Prince, R. A. Mook, R. W. Hendren, and R. W. Wagner, Proc. Natl. Acad. Sci. U.S.A. 
93, 3176-3181 (1996). 

[18] H. Farhood, R. Bottega, R. M. Epand, and L. Huang, Bioch. et Biophys. Acta 1111, 
239-246 (1992). 

[19] D. D. Lasic, H. Strey, M. C. A. Stuart, R. Podgornik, and P. M. Prederik, J. Am. Chem. 
Soc. 119, 832-833 (1997). 

[20] J. N. Israelachvili, D. J. Mitchell, and B. W. Ninham, B. W., J. Chem. Soc, Faraday 
Trans. 72, 1525-1568 (1976). 

[21] P. S. Kuhn, Y. Levin, and M. C. Barbosa, Chem. Phys. Lett. 298, 51-56 (1998). 

[22] P. W. Debye, and E. Hiickel, Phys. Z. 24, 185-206 (1923). 

[23] Y. Levin, Europhys. Lett. 34, 405-410 (1996). 

[24] S. Alexander et. al, J. Chem. Phys. 80, 5776-5781 (1984). 

[25] M. E. Fisher, and Y. Levin, Phys. Rev. Lett. 71, 3826-3829 (1993). 

[26] P. Flory, Principles of Polymer Chemistry (Cornell University Press, Ithaca, New York, 
1971). 

[27] G. S. Manning, J. Chem. Phys. 51, 924-933 (1969). 
[28] G. S. Manning, Q. Rev. Biophys. 11, 179-246 (1978). 

[29] A. V. Gorelov, E. D. Kudryashov, J. -C. Jacquier, J.-C., D. McLoughlin, and K. E. 
Dawson, Physica A 249, 216-225 (1998). 

[30] K. Shirahama, K. Takashima, and N. Takisawa, Bull. Chem. Soc. Jpn. 60, 43-47 (1987). 



8 



FIGURES 
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FIG. 1. Schematic Illustration of the Model: Complexes formed by DNA molecules (rectan- 
gles), counterions (solid circles) and flexible surfactant molecules ( open circles with chains) are 
surrounded by free counterions, colons ( open circles), and surfactant molecules. 
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FIG. 2. Lattice model for forming a cluster. Each monomer can be either occupied by a 
counterion (C) , a surfactant {S), or neither (— ). Interaction between the counterions is purely 
electrostatic, while two neighboring surfactant molecules gain energy due to hydrophobic effects. 




Molar concentration of surfactants 

FIG. 3. Fraction of monomers along the DNA segment associated with surfactant molecules, 
'>^suTf/Z, as a function of bulk density of surfactant, ps- The solid line is derived from the solution 
of Eq.(8), while the circles are the experimental data[29]. The experimental value of the hydropho- 
bicity parameter is x = — 3.5A;bT'[20]. The concentrations of DNA and of added salt are 2 x 10~^M 
and 18mM, respectively. The length of a DNA segment is 220 base pairs. 
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